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Parallel Replica Dynamics 



The system vibrates in 3-N dimensional basin many times before finding an 
escape path. Using many trajectories in parallel, we can find a correct 
escape event more quickly. (I.e., the probability for escaping along a certain 
path should be proportional to the rate constant for that path.) 
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Parallel Replica Dynamics 


Parallelizes time evolution 
Assumptions: 

- infrequent events 

- exponential distribution of first-escape times 



AFV, Phys. Rev. B, 57, R13985 (1998) 
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Parallel Replica Dynamics Procedure 


Replicate entire system on each of M processors. 
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Parallel Replica Dynamics Procedure 


Randomize momenta independently on each processor. 
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Parallel Replica Dynamics Procedure 


Run MD for short time (x dephase ) to dephase the replicas. 
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Parallel Replica Dynamics Procedure 


Start clock and run thermostatted MD on each processor. 
Watch for transition... 
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Parallel Replica Dynamics Procedure 


Stop all trajectories when first transition occurs on any processor. 
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Parallel Replica Dynamics Procedure 


Sum the trajectory times over all M processors. Advance simulation 
clock by this t sum 
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Parallel Replica Dynamics Procedure 


On the processor where a transition occurred, continue trajectory for a 
time x corr to allow correlated dynamical events. 
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Parallel Replica Dynamics Procedure 


Advance simulation clock by x corr . 
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Parallel Replica Dynamics Procedure 


Replicate the new state and begin procedure again. 
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Parallel Replica Dynamics 




The summed time (t sum ) obeys the correct exponential distribution, 
and the system escapes to an appropriate state. 

State-to-state dynamics are thus correct; x corr stage even releases 
the TST assumption [AFV, Phys. Rev. B, 57, R13985 (1998)]. 

Good parallel efficiency if x rxn / M >> x dephase +x corr 
Applicable to any system with exponential first-event statistics 
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Detecting a transition 

best method depends on the system 

simple method for EAM metal systems: 

periodically perform steepest-descent quench; 
see if geometry at basin minimum has changed 



can also use: 

- change in bond connectivity (Kum, Uberuaga) 

- change in local order parameter 






Examples of ParRep studies 



Ag 169 /Cu(100), magic cluster, 
Uche et al, 2009. 



sliding, 500 i^m/s 
Mishin et al, 2007. 




Friction force microscopy, 
Dong et al, 2009, 2010, 
2011. 



Cu void collapse to 
SFT, ^is, Uberuaga 
et al, 2007. 



Hexadecane pyrolysis, 
l^s, Kum et al, 2004. 


Recent brief AMD review: Perez et al, Ann. Rep. Comp. Chem. 5, 79 (2009)i_os Alamos 











































